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[ Abstract 

, We study two-dimensional Af = (2, 2) SU (N) super Yang-Mills theory on Eu- 

clidean two-torus using Sugino's lattice regularization. We perform the Monte-Carlo 
simulation for N — 2,3,4,5 and then extrapolate the result to N = oo. With the 
periodic boundary conditions for the fermions along both circles, we establish the ex- 
istence of a bound state in which scalar fields clump around the origin, in spite of 
' the existence of a classical flat direction. In this phase the global (Zjv) 2 symmetry 

turns out to be broken. We provide a simple explanation for this fact and discuss its 
physical implications. 
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1 Introduction 



Large- N supersymmetric Yang-Mills theories (SYM) are promising candidates for nonper- 
turbative formulations of the superstring theories [IJ EJ [3j HJ [5] . Recent developments of 
discretization technique enable us to study nonperturbative aspects of these theories. (For 
a recent review, see e.g. (6j [7].) Especially, one-dimensional maximally supersymmetric 
gauge theory has been studied extensively [81 [HI [TO], [TTJ H] an d the gauge/gravity dual- 
ity [H [5] has been confirmed very precisely, including the stringy a' corrections [12] . By 
assuming that the gauge/gravity duality holds, then the Monte-Carlo simulation of the 
gauge theory provides a new and powerful tool to study the physics of black holes. 

Next simplest model to study is two-dimensional theory, of which some lattice regular- 
izations are known. Maximally supersymmetric 2d SYM is expected to have dual Dl-brane 
description in type IIB superstring theory [5^1 . By compactifying the spatial direction, one 
obtains Dl-branes winding on the compactified direction. By taking T-dual, one obtains 
a system of DO-branes in compact space, which can have several phases - if DO-branes are 
smeared along the compactified direction it is a black string, and if they are localized in a 
small region it is a black hole. The transition between these phases (the Gregory-Laflamme 
transition [JJ5]) corresponds to the breakdown of the global Zj\r symmetry in the gauge 
theory [II]. By studying the Gregory-Laflamme transition in supergravity and then by 
using the gauge/gravity duality to translate the supergravity to the gauge theory, one can 
study the phase structure of the gauge theory at strong coupling [15] . Using the duality 
in the opposite direction, we can study the detail of the stringy correction to the Gregory- 
Laflamme transition with the Monte-Carlo simulation of the gauge theory. At present, 
it is difficult to study the maximally supersymmetric 2d SYM by using the Monte-Carlo 
simulation. However, four-supercharge system has been studied extensively by using the 
formulations free from fine tuning. 

In this paper, we study two-dimensional M = (2, 2) SU(N) SYM on two-torus using 
Sugino's lattice regularization j!6lpl. In this model, the restoration of supersymmetry with- 
out fine-tuning has been tested~tensively [Eg. We study N = 2, 3, 4, 5 and extrapolate 
to the planar limit N = oo. The action in the continuum is obtained from 4d J\f = 1 SYM 
through the dimensional reduction, and is given by 



where fj, and v run x and y, i and j run 1 and 2, and T 1 = (T^, F l ) are gamma matrices in 
four dimensions. Xj are N x N hermitian matrices, ip a are N x N fermionic matrices with 
a Majorana index a and the covariant derivative is given by Dn = da — i[Au, ■ ]■ The only 
parameters of the model are the size of circles L x and L y . (Note that the coupling constant 
can be absorbed by redefining the fields and coordinates. In other words, the 't Hooft 
coupling A which has mass dimension 2 can be taken to be A = 1. Then the strong coupling 
corresponds to the large volume.) There are several motivations to study this system. 
Firstly, it is the simplest SYM in two dimensions which can be studied nonperturbatively 

1 Of course this model has another interesting interpretation as the "matrix string theory" [3]. 

2 For other regularizations of this theory, see [17] . 

3 See also [191 1201 12T1 122] . For simulations using other formulations, see e.g. [23] . 
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by lattice simulation. Especially, notorious "sign problem" is absent. Secondly, we can 
expect that it is qualitatively similar to maximally supersymmetric (M = (8,8)) SYM, 
which is conjectured to be dual to type II superstring. Thirdly, its bosonic cousin is 
studied in [21] and it is interesting to compare the phase structure. In the bosonic model, 
the Z^r symmetry is broken below the critical volume. In the supersymmetric model it is 
expected to be broken in any finite volume [25]. The argument in [25] is valid only large 
and small volume region and it is desirable to check the breakdown of Zjv at intermediate 
volume. 

An obstacle for the simulation is the existence of the flat direction, along which two 
scalar fields X\ and X2 commute. In contrary to a theory on M 1 ' 3 , there is no superselec- 
tion of the moduli parameter in this case. That is, eigenvalues of scalars are determined 
dynamically. Therefore, some mechanism which restrict eigenvalues to a finite region is 
necessary in order for the stable simulation. One possible way is to introduce an IR regu- 
lator and gradually remove it. In |T8l Ell [22] this method has been applied to the SU(2) 
theory at finite temperature. In those works, a mass term of scalars has been introduced 
and physical quantities are evaluated by an extrapolation to the massless limit or eval- 
uated with small scalar mass. In this case, as we will see, in fact the scalar eigenvalues 
spread as IR regulator is removed. In string terminology, this phase can be understood as 
a gas of freely propagating D-branes. 

At large-iV, there is a more interesting phase, namely a bound state in which eigen- 
values clump to a small region [25J. It is metastable at finite- N and becomes stable at 
large-iV. This bound state is a cousin of the one in one-dimensional system, which has 
been found in [9], and corresponds to the black brane background in type II superstring. 
(The black brane solution in supergravity corresponds to the bound state of the D-branes. 
Because the scalar fields represent the collective coordinates of the D-branes, the bunch 
of the D-branes is nothing but the bound state of scalar eigenvalues.) In this paper we 
study this phase. Because the very existence of this bound state is nontrivial from a 
purely field-theoretical point of view, we provide numerical evidence. As we will see, with 
periodic-periodic (P-P) boundary conditions for fermions, we construct the bound state 
explicitly. In this bound state, there are two possible phases, namely the (Z^) 2 broken 
and unbroken phases. (The model has global (Zjy) 2 symmetry, which shifts the complex 
phases of Wilson loops winding on circles. If Wilson loops are non-zero, then the (Z^) 2 
symmetry is spontaneously broken.) We study the phase structure under the variation of 
the periods. Because of the limitation of the resources, we study only the case of L x = L y . 
We confirm that, in the bound state, the (Z^) 2 symmetry is broken as discussed in [25 1 . 
With antiperiodic-periodic (A-P) boundary conditions, we need rather large N to find 
such an bound state and we could not construct it numerically. Hence we cannot discuss 
the thermal properties of the black brane in this paper. However we can expect that we 
can study the finite temperature system in near future, by using a faster computer. 

This paper is organized as follows. In § [2] we review the conjectured phase diagram 
and explain its physical implications. In §[3] we show the numerical results. We give the 
conclusion and discussion in § HI The simulation details are explained in the appendices. 
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Model 


large volume 


finite volume 


small volume 


4d Af 


= 1, P-P 


(Z^y) 4 and (broken) 


(Z^) 4 and (broken) 


(Zjy) 4 and (broken) 


2d M 


= (2,2), P-P 


{In? 


(broken) 


(broken) 


2d N 


= (8,8), A-P 




(phase transition(s)) 


(broken) 



Table 1: Conjecture for phases of various super Yang-Mills. In this paper, we confirm the 
2d M = (2, 2) case numerically. 



2 Conjectured phase structure 

In this section we review the expectations on the phase structure. In short, the Zjv 
symmetry is broken in the bound state. The phases are summarized in Table [TJ 

2.1 4d N = 1 SYM on Euclidean four-torus 

Let us start by considering the 4d M = 1 SYM on four-torus T 4 . For all circles we impose 
periodic boundary conditions for both bosons and fermions. Let L M (/u = 1, • • • , 4) be the 
periods along four directions and Wn be the Wilson loops winding on each directions, 

W^ = Pexp(iJ^ " dx^A^j , (2) 

where the contraction over fi is not taken in the right hand side. Under the global (Ztv) 4 
transformation, are multiplied by phase factors, 

Wfl _> e ^/N w ^ (n^T). (3) 

Therefore the (1>n) ls broken if TrW^ have nonzero expectation values. In 4d J\f = 1, 
there are (at least) two phases - (Zjv) 4 broken and unbroken. (Note that the Z^v symmetry 
can be broken only in the large- N limit.) 

The existence of the (Ztv) 4 unbroken phase can be shown as follows [26]. Consider the 
situation in which one of the periods L\ is much smaller than other three periods L2,L^ 
and L4. Then essentially we obtain SYM on R 3 x S 1 . If we take L\ <C ^qcd *C L2, L3, L4, 
the effective potential of the Wilson line phase In W\ can be calculated perturbatively and 
it gives exactly zero. By taking monopole and instanton effects |27| into account, it turns 
out that the eigenvalue of the Wilson line (the Wilson line phases) repel each other once 
they spread uniformly. Therefore, Z^-unbroken configuration is stable. This calculation 
itself is done in a specific limit, but due to the supersymmetry the same result should hold 
at any volume [26] . 

The (Zn) 4 broken phase can be found at small volume limit [25j . Suppose that (Z^v) 4 
is broken, or equivalently, the Wilson line phases clump in a small region around the origin. 
If the size of the phase distribution is small enough, then the system can be approximated 
by its zero-dimensional reduction^- This model has been studied extensively. Especially, it 

4 In the dimensional reduction, KK modes are neglected while the effective mass term coming from the 
commutator term TrlA^, A v ] 2 in the field strength term is kept. In order for this approximation to make 
sense, the effective mass must be sufficiently smaller than the KK mass, that is, eigenvalues of must 
be small enough. Then the Wilson loop is close to 1 and hence the Zjv symmetry is broken. 
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is known to have a bound state of eigenvalues even though it has a flat direction classically, 
and the size of the eigenvalue distribution is known EH] . It is small enough so that 

the assumption that the system is described by the zero-dimensional reduction is correct. 
Therefore, a (Z^r) 4 broken phase exists at small volume. We expect this phase persists at 
any finite volume due to the supersymmetry. 

2.2 2d TV = (2, 2) SYM on two-torus 

2d J\f = (2, 2) SYM is obtained from 4d N = 1 through the dimensional reduction. When 
we take two circles in 4d theory to be small, in order for the dimensional reduction to 
work the Zjy symmetries along these circles must be broken. In such a phase Wilson line 
phases clump, so it corresponds to the bound state in 2d theory. Therefore, it is natural 
to assume that the bound state in 2d theory corresponds to (Zjv) broken phase in 4d 
theory. If this is the case, the (Z^r) 2 symmetry of the 2d theory should be broken. 

At large volume, (1>n) 2 symmetry should be restored. This is because (1>n) 2 is essen- 
tially continuous U(l) 2 symmetry at large- iV and hence at noncompact two-dimensional 
space it is unbroken according to the Coleman's theorem. 

In § O we will confirm the above statement by numerical simulation. The discussion 
above does not exclude an existence of a Z^v unbroken phase; see a remark in § 13.11 

2.3 2d W = (8,8) SYM : Z N broken phase as a black hole 

For 2d N = (8, 8) SYM, there is a dual gravity interpretation. Here we consider the 
correspondence between thermal SYM and black branes in type II supergravity. First, let 
us briefly describe the SYM in terms of D-branes. We impose the antiperiodic boundary 
condition, for temporal direction (we take y to be temporal direction). In this case, the 
bound state (i.e. a state in which scalar eigenvalues clump around the origin) is dual to 
the system of coincident N Dl-branes at finite temperature. The temperature T is inverse 
of the (3 = L y . By taking T-duality along the spatial circle (x-direction), we obtain the 
system of N DO-branes, sitting in the compactified spatial dimension of radius (2ir) 2 /L x . 
The Wilson line phases correspond to the positions of the DO-branes along x-direction. 

According to the gauge/gravity duality conjecture [5], at low temperature this system 
is well approximated by type II supergravity, and for fixed (3 there is a first-order phase 
transition associated with the breakdown of the Z n symmetry along x-direction [14|, [T5] . 
Below the critical value of the radius L XjC ~ yj]3, Zjv is broken. Along the temporal 
direction, Z^v is conjectured to be broken at any nonzero temperature, that is, the system 
is always deconfined. For Id system obtained by dimensional reduction from 2d J\f = (8, 8) 
SYM it has been confirmed by using the Monte-Carlo simulation ffi 112} ITT]. 

The physical interpretation of the Z^r breakdown in string theory is simple. In Z^r 
unbroken phase, DO-branes (Wilson line phases) fill the compact direction uniformly. It 
is a "uniform black string"^. If Zjy is broken completely, DO-branes clump to a small 
region along compact direction; it is similar to the usual "black hole" (black 0-brane). 
This black hole is small at low temperature, and cannot wind on the spatial direction. At 

5 Do not confuse with "black 1-brane" which is a classical solution to type IIB supergravity corresponding 
to a bunch of Dl-branes, while the black string is a solution to type IIA. The black string and the black 
1-brane are related by the T-duality. 
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high temperature the black hole becomes large, and gradually fills the compact dimension, 
wraps on it and becomes the black string. At intermediate temperature the "nonuniform 
black string" may exist. In this phase, "Lm is broken and phases of the Wilson line is 
distributed nonuniformly, but still the density is nonzero everywhere. 

In the high-temperature limit, fermions decouple and the system reduces to Id bosonic 
Yang-Mills. In [31] this limit was studied by the Monte-Carlo simulation and it was found 
that the nonuniform black string phase exists indeed. By fixing the temperature and 
varying the size of the spatial circle, there are two phase transitions, namely between 
uniform and nonuniform black strings and between the nonuniform black string and the 
black hole. The orders of the transitions are of third and second order, respectively. 

The situation is similar for (p + l)-dimensional SYM with 16 supercharges, in which 
the bound state is dual to a bunch of Dp-branes. At low temperature we can analyze 
it using type II supergravity and (Ztv) p symmetry broken/unbroken phases correspond 
to smeared DO-brane solution and black 0-brane solution. For detailed phase structure, 
see [52]. 

We emphasize again that the bound state of scalar eigenvalues is necessary for the 
black brane description. In order to understand this statement, let us remind that the 
black brane is the bound state of very large number of D-branes (indeed the gravity picture 
is valid in the large- N limit, where A is nothing but the number of D-branes). Because 
the scalar fields represent the collective coordinates of the D-branes, the bound state of 
the D-branes is exactly the bound state of scalar eigenvalues. 

2.4 Zjv breakdown and volume (in) dependence 

Before concluding this section, we argue general features of the large-A volume indepen- 
dence. 

In the large-A limit, if the Z^v symmetry is not broken then the system is volume 
independent (Eguchi-Kawai equivalence) [331 123] • That is, physical quantities do not 
depend on the volume up to a trivial proportionality factor. For example, the free energy 
per unit volume is volume- independent. This property is practically very useful when one 
studies the large-A field theories numerically, because by using small-volume (sometimes 
zero- volume) lattice one can save computational cost. Recently its application for a study 
of the Ztv unbroken phase of 4d M = 1 pure SYM [26] has been discussed [331 [35]. In the 
case of 2d N = (2,2), however, the Z^v symmetry is broken (although we cannot prove 
the non-existence of the Z^v unbroken phase) and hence the theory is volume dependent. 

There is another formulation utilizing the Eguchi-Kawai equivalence [36J. In this con- 
struction, matrix quantum mechanics around a certain background is equivalent to 4d 
gauge theory on M x S 3 . This theory is manifestly volume-dependent because the curva- 
ture of the sphere emerges as a parameter. This technique can be used to formulate 4d 
J\f = 1 SYM [37] • Similarly, 3d gauge theory on Rx5 2 can be formulated by expanding the 
matrix quantum mechanics around a fuzzy sphere and then taking the commutative limit. 
We can also formulate 3d theory on S 3 and 2d theory on S 2 by using a zero-dimensional 
matrix model. These models are presumably corresponding to the Z^v broken phase. 
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3 Numerical results and the phase structure 



In this section we study the phase structure of 2d N = (2, 2) SYM in the case of P-P 
boundary condition numerically. Because of the limitation of the resources, we study only 



the case that L x = L y = L. 



3.1 The existence of the bound state 



As we have explained, the system should have a bound state in which scalar eigenvalues 
clump around the origin. Although this is stable at large-A^, for small N it is at most 
metastable; it often collapses and eigenvalues spread along the flat direction. In order for 
the stable simulation, we add the mass term for scalars to regularize the flat direction, 

S = N [ dx [ dy Tr (-F 2 + -(D^X 1 ) 2 - -[X\X^] 2 + m 2 (X i ) 2 V(fermionic part), 
Jo Jo \4 2 4 / 

(4) 

and gradually remove it. Furthermore we impose the periodic-periodic boundary con- 
dition, because the metastable state becomes more stable since the fermion zero-modes 

provide an attractive force between eigenvalues |38j. In Fig. [H the mass-dependence of 

l 



jjj^j f dx f dy Tr Xfjj is plotted. (The physical volume is L = 0.707.) As we can 
see, it converges to a nonzero value, which suggests the existence of the bound state. As we 
will show in § 13.21 the mass-dependence of the Wilson loops disappears at m < 0.4. This 
suggests that the mass term is so small that it is effective only when the eigenvalue of the 
scalar fields deviates significantly from the typical size of the bound state. Furthermore, 
as we will show in § 13.31 this state corresponds to the bound state in the zero-dimensional 
matrix model. Therefore we conclude that we have constructed the bound state mentioned 
in the previous section. We also plot the same quantity with A-P boundary condition. It 
is clear that the scalar eigenvalues diverges in the m — > limit. Hence this phase is not 
the bound state. Note that the unbounded state exists also with P-P boundary condition. 
Indeed, we observed that the norm of the scalar fields sometimes blows up if we set m = 0. 
This can be interpreted that a meta-stable bound state decayed into an unbounded state. 
For A-P boundary condition, we expect that larger value of N is needed to stabilize the 
bound state, as in Id theory studied in [9]. 

A few remarks on the previous numerical simulations for SU (2) gauge group [181 EH 
[20l [2T| [22] are in order here. In |18^ I2T] [22] the same regularization by introduction of 
the mass term is used. In these works the scalars seems to diverge in m — ► limit (it is 
the same as the one with A-P boundary condition, shown in Fig. [TJ See also |39|.) and 
hence it is plausible that the extrapolation to m = picks up the different phase from 
the bound stat^|. It is plausible that this phase corresponds to the Zjv-unbroken phase 
discussed in [261 ES]. On the other hand, the phase studied in [HI [20] should be the bound 
state, because in these papers the configurations have been made by using the bosonic 
part of the action and the effect of the fermionic part has been taken into account by 
the reweighting method. (Note that in the bosonic model the flat direction is lifted by 
quantum corrections and hence only bound states exist.) 



6 However, the configurations may have a large overlap with the bound state as well since the scalars 
are gathered around the origin by the effect of the finite mass. 
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1.2 




° P-P, SU(2) i — i — i 

0.2 - P-P, SU(3) : ----x ---: - 

P-P, SU(4) i * ! 



0.2 0.4 0.6 0.8 1 1.2 

m 

Figure 1: The mass-dependence of scalar eigenvalues at L = 0.707 is plotted. The vertical 

axis represents ((yj^z f dx J dy TrXf\\ , which is the inverse of the typical value of 
the scalar eigenvalues. In P-P boundary condition the scalar eigenvalues remain finite in 
the m — > limit, while it diverges in A-P boundary condition. Fitting lines are quadratic 
polynomials of m. 



3.2 Large- A" behavior of the Wilson loop and breakdown of Z N symme- 
try 

In this subsection, we study the Wilson loop Tr Wa/N, where Wn is a unitary matrix 
which is defined by 

W x (y) = Pexp fi J dxA x (x,y)^j , W y (x) = Pexp (i J dyA y (x,y)^j . (5) 

On lattice, it can be obtained by multiplying the link variables. Because of the transla- 
tional invariance of the model, the expectation value does not depend on the coordinate. 

At finite N, there is a tunneling between different vacua related by a multiplication 
of a phase factor ([3j) and hence the expectation value (-^TrW^) is zero. (In the large- A 
limit, the tunneling is suppressed and the expectation value can be non-zero.) Therefore, 
in order to see a possible symmetry breakdown at large- A, we fix the (Z^r) 2 symmetry 
so that ReJ2j=iTrW$' becomes maximum. That is, Wjf' is replaced by e 2mn ^/ N w^ 
with a suitable which satisfies the above condition. In our numerical simulation, we 
have performed this Z^r fixing at each measurement. 

Below we show the expectation values of the Wilson loops. We found that (-^ Tr W x ) 
and (jjTrWy) are the same in the error. Therefore we use the average of them (±TrW), 
where W = (W x + W y )/2. (In the bosonic models, they take different values for some 
parameters [2H 132] .) 

In Fig. [21 the mass-dependence of the Wilson loop is plotted. (The size of the torus 
is L = 0.707.) The mass-dependence disappears at small values of m, which allows us to 
use small fixed value of m. In practice, we use m = 0.2, 0.3 and 0.40 

7 Without introducing the mass m, scalar eigenvalues stays around the origin for a while and then go 



8 




Figure 2: The expectation value of the Wilson loop (Tr W/N) at L = 0.707 is plotted for 
< m < 1. At small values of m, the mass-dependence disappears. 



3 0.4 




Figure 3: The expectation value of the Wilson loop {Tr W/N) for N = 3, 4 and 5. Mass 
parameter is m = 0.2. 



0.8 
0.6 
0.4 
0.2 



0.588/N +0.845 
-2.13/N^+0.695 



0.02 0.04 0.06 0.08 0.1 0.12 

1/N 2 



Figure 4: iV-dependence of the coefficients in (|6|), N = 3,4 and 5. Mass parameter is 
m = 0.2. 
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In Fig. [3] we plot the expectation value of the Wilson loop against the size of the torus 
L. For each N, the expectation value can be fitted by 

(^l-TrW^j = c x exp(-d-L), (6) 

where c and d are real and positive. By using the data at L = 0.471, 0.707, 1.061, 1.414 and 
m = 0.2, we obtain c = 0.911 ± 0.012, d = 0.464 ± 0.018 for N = 3, c = 0.877 ± 0.015, d = 
0.538 ±0.022 for N = 4 and c = 0.872 ± 0.016, d = 0.624 ±0.020 for N = 5. Extrapolating 
these coefficients to N = oo, we obtain c = 0.845 ± 0.009 and d = 0.695 ± 0.040 (Fig. SJ). 
Therefore, in the large-iV limit, the expectation value of the Wilson loop is nonzero at 
finite volume, and hence the "Lm symmetry is broken. Note that the TL^ symmetry is 
restored in the large volume limit as expected. 



3.3 Zero-volume limit 

As we have seen above, the simulation data suggests that the Ztv symmetry is broken 
completely in the zero-volume limit, that is, the expectation value of the Wilson loop 

TrW^j becomes close to 1. Then the model reduces to the one-point reduction of 4d 
which has been studied numerically in [29]. Hence, by comparing the 2d 
and Od models, we can check the validity of our simulation. 

Let us start with the action in two dimensions, 

S2d = N £ dx £ dy Tr Q(Z^) 2 - \[Y i: Yj] 2 + -^F 2 ^ . (7) 

Here we denoted the scalar fields as Yi in order to distinguish them from the scalars Xj 
in the zero-dimensional model. (For simplicity, we write down only the bosonic part.) 
By neglecting nonzero- modes, this action reduces to the zero-dimensional one. By using 
X I = (X fl ,X i ) defined by 

Xi = L~^ 2 [ dx I dyY, X il = L-*l 2 [ dx f dy (8) 
Jo Jo Jo Jo 

we obtain the zero-dimensional model with the canonical normalization, 

So^-jTrlXj^j] 2 . (9) 
Here, the gauge field can be obtained by 

\[-i{W^-l)+c.c]~ f dx^A^, (10) 
1 Jo 

to infinity. It is consistent with the interpretation that the bound state is metastable. It is not impossible 
to evaluate expectation values by using only configurations from metastable configuration, although the 
error is rather large and the result is less reliable. The result is consistent with the constant fitting shown 
in Fig. E] 

8 Unless the Zjv symmetry breaks completely, the zero-volume limit is different from the naive dimen- 
sionally reduced model; see [13] . 
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where in the r.h.s. the index /i is not contracted. When we consider the k x k square 

(1) (k) 

lattice, we can calculate k Wilson loops along both x and y directions, , • • • , . 
can be obtained aa^| 



There is a subtlety in the discussion above - "zero-mode" is not a gauge-invariant 
notion. (A constant field configuration, which corresponds to the zero-dimensional model, 
can be transformed to rapidly varying configuration just by a gauge transformation!) 
Therefore, we have to choose a suitable gauge in which zero modes reproduce the zero- 
dimensional matrix model [10]. In the current setup, we should maximize (resp. minimize) 
the zero-mode (resp. nonzero-mode) contributions, so that the configuration is as static 
as possible. For that purpose we choose the gauge so that Tr(Xf + X 2 ) is maximum. 

In this gauge we can compare the small volume behavior with zero-dimensional matrix 
model. The simplest quantity is ( — Tr[Xj, Xj] 2 }, which can be evaluated exactly [41J: 

~Tr[X u X 2 ?) = ^Tr[X 3 ,X^ = \ (l - i,) . (12) 

In Fig. [5] and Fig. EJ the corresponding quantity in two-dimensional SYM is plotted. In 
Fig. [7J and Fig. [SJ / \J Y^i=i Tv" ^ r ~^tft) an< ^ ^he corresponding quantity in two dimensions 
are plotted. In these plots, the results at m = is obtained by extrapolating results at 
m = 0.2, 0.3 and m = 0.4 by a straight line. Then we have fitted them by a polynomial 

A+BL+CL 2 . We have evaluated ^ \JY^i=i ~k Tr X 2 ^ numerically in Od theory and found 
that it is 1.78 ± 0.03 for SU(3). In Fig. (SJ Fig. [6] and Fig. M the data is consistent with 
Od model results for SU(3). (In Fig. [7J the extrapolated value differs slightly from the Od 
result. One possible reason is this quantity is sensitive to the error in the approximation in 
(llip . which is exact in the limit Wn = In-) Since our 2d simulation is smoothly connected 
to the Od model, we conclude that the bound state we have constructed is exactly the one 
discussed in 



3.4 Supersymmetry 

It is important to study whether the bound state preserves supersymmetry or not. Here 
we show the expectation value of the action. In the Sugino model the action is of the form 

S = {Q,0}, (13) 

where Q is one of four supercharges which is exactly kept in the regularization. Therefore, 
the expectation value of the action must be zero if the vacuum is invariant under the 

9 Because X M should be traceless, we have to fix the (Zjv) 2 symmetry so that Re Ysj=i Wf? becomes 
maximum, or equivalently the argument of ^2j =1 TrWff is minimum. (It is the same as the Zat fixing 
performed in § 13.21 ) Equations (|10[1 and hold only with this condition. In our numerical simulation, 
we have performed this Zjv fixing at each measurement. 
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Figure 9: Mass-dependence of the action at L = 0.707. The points at m = are obtained 
by extrapolations. In the m — > limit, where the continuum theory becomes super- 
symmetric, the expectation value is consistent with zero ( SU(4) case seems nonzero but 
consistent with zero in 2 standard deviations). 



supersymmetry generated by Q. As we can see from Fig. [9l the expectation value of 
the action is consistent with zero. 

Strictly speaking, there is a subtlety for studying the breakdown of the supersymmetry 
with periodic boundary conditions. The reason is as follows [20]. In the simulation, we 
obtain the expectation values normalized as 

(o\{Q,Q}\o) , u , 
(0|0) 1 ) 

where the denominator (0|0) is the partition function. In order for the simulation to make 
sense, the denominator must be nonzero, but this condition can be broken. Indeed, if the 
continuum spectrum is absent and the Witten index is well-defined, then the partition 
function with the periodic boundary conditions is nothing but the Witten index, which 
is zero if the supersymmetry is spontaneously broken. In the present case, because the 
continuum spectrum exists due to the existence of the unbounded state and there is an 
ambiguity for the definition of the Witten index, the above argument cannot be applied 
straightforwardly. Although it is difficult to exclude the possibility that the partition 
function is zero, we believe the partition function is nonzero because this system does not 
suffer from the sign problem. Clarification of this point is desirable. 

In order to see whether the supersymmetry is spontaneously broken or not, the simplest 
and unambiguous way is to put the theory at finite temperature and calculate the energy 
density [20J. Supersymmetry is not broken if and only if the energy density is zero at 
zero temperature. It is an important future problem to study 2d theory at large-JV and 
at finite temperature, and confirm that the energy converges to zer 

10 Because we are picking up fluctuations around one specific state, Q-exact quantities can have nonzero 
expectation values without introducing any external fields nor temperature. Remember that we are study- 
ing the bound state only while there is an unbounded state as well. 

11 For matrix quantum mechanics which is obtained from 2d Af — (2, 2) SYM through the dimensional 
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gauge group 


volume 


lattice size 


SU(2) 


0.707 


4x4 


SU(2) 


1.414 


8x8, 6x6, 5x5, 4x4 


SU(3) 


0.354, 0.471, 0.707 


4x4 


SU(3) 


1.061, 1.414 


6x6 


SU(4) 


0.354, 0.471, 0.707 


4x4 


SU(4) 


1.061, 1.414 


6x6 


SU(5) 


0.354, 0.471, 0.707, 1.061, 1.414 


4x4 



Table 2: A list of lattice sizes used in this work. 



As we will see in § 13.51 it is plausible that the explicit super symmetry breaking lattice 
artifacts disappears in the continuum limit. Then, the fact that the expectation value of 
the action is zero suggests the absence of the spontaneous supersymmetry breaking. It is 
desirable to check it more rigorously by measuring the energy of the system. 



3.5 Convergence to the continuum limit 

In this subsection, we show that the lattice spacing used in the current simulation is small 
enough to study the system quantitatively. 

The lattice sizes used in this work are listed in Table. In Fig. [TUl Fig. [TT] and 
Fig. H2J we plot the lattice spacing dependences of the extent of the scalar eigenvalues, 
the Wilson loop and the action. Lattice size was taken to be 4 x 4, 5x5, 6x6 and 
8x8 and other parameters are taken to be N = 2,L = 1.414, m = 0.2. It turns out that 
the expectation values, especially that of the Wilson loop, are not sensitive to the lattice 
spacing used in the simulation. Note that 4x4 lattice at L = 1.414 is the coarsest one in 
this work. Therefore the numerical data used in the previous sections is sufficiently close 
to the continuum limit. 

In order to see that the supersymmetries which are broken by a lattice artifact are 
restored in the continuum limit, we utilize the discrete symmetry of the system which 
is related to the R-symmetry. The action in the continuum has the following discrete 
symmetries (we follow the notation in the Appendix): 

X^-|, ipo^^l, ipi -> -ipo (15) 

and 

2^V>i, X^V'o, ipi -> 2> tPo^X, $->-(/>. (16) 

Because of these symmetries, in the continuum limit, all four Yukawa interaction terms 
Cpi, £f4 (fo r explicit forms, see Appendix) should give an identical expectation value. 
In Fig. [T3l we plot them for the coarsest (L = 1.414 with 4x4 lattice) and the finest 

reduction, Smilga conjectured that the supersymmetry is broken in the bound state phase |42j . In the 
unbounded state, the supersymmetry is argued to be unbroken. The corresponding phase in 2d theory has 
been studied in |21l 1181 12"2"] and the ground state energy was found to be consistent with zero |22j . 
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Figure 11: Dependence of the Wilson loop to the lattice spacing at N = 2, L = 1.414, 
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Figure 13: Expectation values of Lpi{i = 1, ...,4), the almost degenerate behavior implies 
the that the simulation is close to the continuum limit. The lowest plot uses the same 
configuration used to plot Fig. 8 in 



(L = 0.354 with 4x4 lattice) cases with N = 5. El We also plot N = 2 case using the same 
configuration used to plot Fig. 8 in [18]. The plots show that the expectation values from 
the finest lattice are almost degenerated for all Cfi {i = 1,-4) and it strongly suggests 
that this discrete symmetry is restored in the continuum limit. Therefore we can expect 
that supersymmetries broken by a lattice artifact are restored in the continuum limit, 
because they are related to the exactly kept supersymmetry via this discrete symmetry. 
Note that, although the coarsest case would have some effects from the lattice artifacts, 
at least for the quantities like the Wilson loop such effects are negligible, as we have seen 
above. 

A remark on the previous simulation [18] is in order here. In [18], the restoration 
of supersymmetries which are explicitly broken by lattice artifacts has been studied. In 
that work, expectation values of two-point functions are used. (In the present simulation, 
this technique cannot be used because the lattice is too small to calculate the two-point 
functions.) In Fig. 8 and 9 in [TS], four different two-point functions are plotted for P-P 
case which should be degenerated because of the above discrete symmetries. The plots, 
however, does not show the degenerate behavior and hence it seems that the simulation is 
far from the continuum limit The non-degenerate behavior is again found in the lowest 
plot in Fig. [13] while in the first plot from the current simulation almost degenerated 
behavior can be seen. 



4 Conclusion and discussions 

In this paper we have studied large- N properties of two-dimensional N = (2, 2) SYM. 
Especially we have established the existence of the bound state in which scalar eigenvalues 

12 Because we use a unit 't Hooft coupling A = Ng 2 — 1, and because of a difference of the normalization 
of the kinetic term, the lattice spacing defined by the same manner as [18] becomes smaller by a factor 

y/Nj2. 

13 Because of this and noisy results for partially conserved supercurrent relation, [18] did not study P-P 
case extensively. 
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clump around the origin. It makes the simulation well-defined in spite of the existence of 
the flat direction along which scalar eigenvalues spread to infinity. We also have shown 
numerically that at finite volume global (Z^r) 2 is broken. This symmetry is restored in 
the large volume limit. 

If the (Z^r) 2 symmetry were not broken, then because of the Eguchi-Kawai equiva- 
lence the expectation values of the loops would not depend on the volume in the large- N 
limit \33\ I24j . The Eguchi-Kawai equivalence works for bosonic Yang-Mills only above 
a critical volume [23] because the (Z^) 4 symmetry is broken below it. The twisted [33] 
and quenched [33] Eguchi-Kawai models were believed to cure this problem, but recently 
turned out to fail in the large- N limit [451 |4"6] . The reason is that the backgrounds col- 
lapses due to large fluctuations. (Another deformation has been proposed in |47|.) For 
supersymmetric theories the Z^v symmetry was expected not to be broken, but as shown in 
this paper, it is not necessarily the case. However, once we combine the idea of twisted or 
quenched Eguchi-Kawai model with the supersymmetry then the Z j\r-unbroken background 
becomes stable 148 1 and the Eguchi-Kawai equivalence can hold. A concrete proposal has 
been given in [36] lj. Studying this direction further is very important, because using the 
Eguchi-Kawai equivalence we can regularize planar SYM in higher dimensions (3d and 4d) 
without using lattice, and hence we may avoid the difficulties in lattice SUSY. 

Although finite-temperature properties of the bound state is very interesting in con- 
nection to the black hole thermodynamics, it is much more difficult to study because we 
need to take N to be rather large. The same difficulty exists in matrix quantum me- 
chanics, though large enough N (say N = 16) can be taken in that case [9] because the 
model requires less computational resources. In [llj . to avoid this difficulty at rather 
small N, finite-temperature properties of the matrix quantum mechanics have been stud- 
ied by performing the simulation with the periodic boundary condition and then taking 
into account the effect of the antiperiodic boundary condition by the reweighting method. 
This method would work for the Z^y-broken phase ("black hole" and "non- uniform black 
string"). However, for the Zjv-unbroken phase, this method might not work because we 
expect a severer overlapping problem. In any case, we expect an unambiguous study of 
the thermal properties with large enough N will be possible in near future. It will provide 
us with valuable insights into black hole/black string physics. 
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A Simulation details 
A.l The Sugino model 

We consider the N = (2,2) supersymmetric Yang-Mills theory on T 2 , whose action is 
given by 

^continuum 

= N £ X dx£ V dy Trl-F^ + ^D^) 2 - ^,X,f - ^VDrf - -^PpQ,^;]} , 

(17) 

which is obtained from four-dimensional M = 1 SYM through the dimensional reduction. 
As a discretization, we use Sugino's lattice action [16] l 15 l . 



f 3 - 6 - \ 

^lattice — Qx&y y~] < y^CBijx) + y^ y CFi(x) > + (auxiliary field), 

x U=l i=l J 



where 



a%a y 



x ((f>(x) - U x (x)(f>(x + a x x)U x {x) l ) } 



+ Tr{(<j)(x) - U y (x)<f)(x + a y y)U y {x)- 1 ) 

CL X (ly 



and 



(18) 



C B1 (x) = -^Tr[0(x),0(x)] 2 , (19) 
N n 

C B 2(x) = -^Tr<f> TL (x) 2 , (20) 
-CbsO?) = -3— Tr{(^(x) - U x {x)4>(x + a x x)U x (x)~ l ) 



x (<£(x) - U y {x)4(x + a y y)Uy{xY l ) } (21) 



iV 

£fi x = -jnTr r? x ,77* , (22) 

Aa x a y 

N 

C F 2{x) = -—Tr{ X {x)[4>{x),x(x)]), (23) 
a x a y 

Cfs(x) = --^Tr{Mx)Mx) (Hx) + U x (x)Hx + a x x)U x (x)- 1 )} , (24) 
a x a y 

C F4 (x) = ^3 Tr {Vi(x)Vi(x) (0(x) + Uy(x)0(x + a y y)Uy(x)- 1 )} , (25) 

Ox 0^ 



1 Here we follow the notation in [191 118] . Under a suitable representation of the Gamma matrices, the 
fermion ip can be taken as ip T = (^o, V'l) X: 
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£fs(x) = i-j-2 Tr (x(x) ■ Q6(x)J , (26) 

£f6(£) = 2Y{-0o (r?0?) - U x (x)r](x + a x x)U x (x)~ l ) } 

-»— 3 Tr{^i (v(x) - U y {x) V (x + a y y)U y (x)- 1 ) }, (27) 

where U(x,/i) are gauge link variables, 4>{x) is a complex scalar, r/(x), x(af) and ip^x) are 
fermion field, and a y are lattice spacingsEl, e is a real parameter which must be chosen 
appropriately for each TV" (in this work, we used e = 2.6), 

= T-T-n^fv^ ' * tl{s) = " ^ ( Tr ' (28) 

where P(x) = U x (x)U y (x + x)U x {x + y)Uy(x) is the plaquette variable, and Q generates 
one of the four supersymmetries, 



QU^x) 


= iip^U^x), 


(29) 




= iip^(x)ip^(x) - i(<j)(x) - U^(x)4>(x + a^U^x)" 1 ), 


(30) 


QHx) 


= o, 


(31) 


Qx(x) 


= H(x), 


(32) 


QH{S) 


= [4>(x),x(x)], 


(33) 


Q4>(x) 


= r)(x), 


(34) 


Qrj(x) 


= [4>(x), </>(£)]■ 


(35) 



Sugino's action Siattice is invariant under the supersymmetry generated by Q, because Q 
is nilpotent up to gauge transformation and S can be written in a Q-exact form. 

In [16] . using super-renormalizability and symmetry argument, it was shown that other 
three supersymmetries, which is broken by a lattice artifact at the discretized level, is 
restored in the continuum limit. Furthermore, in [18] . this restoration has been confirmed 
explicitly by the Monte- Carlo simulation. 



A. 2 Simulation 

We have adopted the rational hybrid Monte-Carlo algorithm [50J. We have use the 
code [51] based on the Remez algorithm to find necessary coefficients in the simulation. 
In 2d J\f = (2, 2) SYM, the complex phase of the fermion determinant is absent in the 
continuum limit and at discretized level only small phase appears as a lattice artifact. 
In this work we have ignored it. Fermi QCD/MDP [52] has been used to develop the 
simulation code. 

Because of the limitation of the resources, we have concentrated on the square torus, 
L x = L y = L. We took the number of sites and lattice spacings in two directions to 
be the same. For each set of parameters, we have collected 1000 — 2000 samples of 
configurations. We have evaluated the error by using the Jack Knife method and it turned 
out the autocorrelations are sufficiently small. 

16 In the actual simulation we have used the isotropic lattice, a x = a v . 
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